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g : ABSTRACT 
(N 

We present ultraviolet and optical spectra of the first bright spot (PA = 29°) on 
^ . Supernova 1987A's equatorial circumstellar ring taken with the Space Telescope Imaging 

^ ! Spectrograph. We interpret this spot as the emission produced by radiative shocks that 

occur where the supernova blast wave strikes an inward protrusion of the ring. The 
^ _ observed line widths and intensity ratios indicate the presence of radiative shocks with 

' velocities ranging from 100 to 250 km entering dense 10^ cm~^) gas. These 

■ observations, and future observations of the development of the spectra and line profiles, 

■ provide a unique opportunity to study the hydrodynamics of radiative shocks. 
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1. INTRODUCTION 

After the discovery of the circumstellar ring around Supernova (SN) 1987A, several authors 
pointed out that the supernova blast wave should strike the ring some 10 - 20 years after the 
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supernova explosion (Luo & McCray 1991; Luo, McCray, & Slavin 1994; Chevalier & Dwarkadas 
1995; Borkowski, Blondin, & McCray 1997b). The estimated time of first contact was uncertain, 
mainly because it depended on the unknown density of circumstellar gas between the supernova 
and the ring. 

This event is now underway. In 1995 an unresolved spot at position angle ~ 29° (Spot 1) on 
the near (North) side of the ring began to brighten (Garnavich, Kirshner, & Challis 1997). This 
spot continues to brighten, and within the last year several new spots have appeared (Lawrence & 
Crotts 2000; Garnavich, Kirshner, & Challis 2000; Lawrence et al. 2000). A shock interpretation 
for these spots was suggested by Space Telescope Imaging Spectrograph (STIS) observations of 
the ring taken in 1997 (Sonneborn et al. 1998). The spectral image of the ring in Ha showed a 
streak at the position of Spot 1 indicating blueshifted emission with velocities up to ~ 250 km s^^. 
Evidently, at the positions of the spots, the blast wave is striking inward protrusions of the dense 
equatorial ring. 

Here (§2) we report on early optical and ultraviolet STIS spectra of SNR 1987A which show 
several broad lines associated with Spot 1. In §3 we present the radiative shock model for the 
emission from the spots, and in §4 we discuss what we can infer from the spectrum and its evolution. 

2. OBSERVATIONS 

The Supernova INtensive Study (SINS) collaboration obtained STIS spectra of Spot 1 at 
ultraviolet and optical wavelengths in 1997 September and 1998 March, respectively. The ultraviolet 
spectrum (not shown) was taken on September 27, 1997 (3869.3 days since explosion) with the 
G140L low dispersion grating (1150 - 1736 A, Av 300 km s~^) and a 0'.'5 slit. Five spectral 
images totaling 11,200 s were combined. Emission lines from nearly stationary gas in the ring are 
apparent at Nv AA1239, 1243, Oiv] A1400, Niv] A1487, and Hell A1640. Emission from Spot 1 is 
also visible in these lines, as well as at Siiv AA1394, 1403, [Niv] A1483, and Civ AA1548, 1551. The 
low- velocity Si IV and C IV emission from the ring is blocked by interstellar absorption (Pransson et 
al. 1989). Due to geocoronal Lya emission no information about Lya from the spot was obtained. 

The optical spectrum was taken on March 7, 1998 (4030.0 days since explosion), with a 0'.'2 slit 
and the G750M grating (6295 - 6867 A, Av « 50 km s"'''). Three spectral images totaling 8,056 s 
were combined. A section of the optical spectrum is shown in Fig. 1. The central horizontal streak 
(which extends to ±3, 000 km s~^) is Ha emission from the radioactively heated supernova 
debris. Pairs of bright spots where the slit intersects the stationary inner ring are evident in [N li] 
AA6548,6584 and Ha (not displayed but also present are [Oi] AA6300,6364, He I A6678, and [Sii] 
AA6717, 6731). Three fainter spots due to the outer rings are also visible in Ha and [Nil]. Broad 
streaks of emission from Spot 1 appear on the upper segments of the ring lines. These streaks are 
inside the ring emission indicating that Spot 1 lies at the inner edge of the ring. 
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Fig. 1. — (a) Slit orientation on image of SN 1987A's triple ring system, (b) Section of STIS G750M 
optical spectrum (the intensity scale has been stretched to show the hot spot emission). 



The spectra of Spot 1 are extracted from the STIS spectral images. For the optical emission 
lines, the major source of contamination is the emission from the nearly stationary [v = 10.5 it 
0.3 km s"\ Crotts & Heathcote 1999) circumstellar ring filling the Qf.'2{^ 100 km s~^) slit. The Ha 
line profile is shown in Figure 2. The parts of the line profiles that lie outside of the ring emission 
are fit well by Gaussians. The velocity shifts of the Gaussian fits are poorly constrained since the 
spot's spatial position in the slit is uncertain. In the ultraviolet the ring emission is much weaker 
and therefore its contribution to the extracted Spot 1 spectrum is neglected. The UV lines are 
poorly resolved at the resolution of the G140L grating. It is possible to account for the detector 
line spread function in order to obtain intrinsic widths for the UV lines (Pun et al. 2000). While 
this method introduces errors in the determined width, we find that all the UV lines have intrinsic 
widths < 400 km s~^. 




Fig. 2. — Ha profile of Spot 1 extracted from the G750M spectrum. Low velocity emission (shaded 
region) from unshocked stationary ring material defines the standard of rest. A Gaussian fit to the 
broad emission from shocked gas is shown. 



All the lines from Spot 1 were fit with Gaussians in order to extract their fluxes and full widths 
at half maximum (FWHM) (see Table 1). The errors stated in Table 1 are la statistical errors 
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to the fits. Close multiplets (Nv AA1239, 1243 and CiV AA1548, 1551) were fit simultaneously by 
requiring a set wavelength separation, identical widths, and line ratios dictated by their oscillator 
strengths. For the Nv doublet, the emission from Spot 1 overlapped with very broad Lya emission 
(~ ±15,000 km s"^) from the reverse shock (Sonneborn et al. 1998; Michael et al. 1998). This 
Lya background was removed using a quadratic fit and is the source of the large uncertainty in 
the Nv AA1239, 1243 flux. At the G140L resolution, the Si IV A1403 emission of Spot 1 is mixed 
with that from the Oiv] A1400 multiplet. Using the measured Si IV A1394 flux and assuming 
/(1394)//(1403) = 2, the Si IV A1403 flux was determined and removed from the Oiv] A1400 
multiplet. 



Emission Lines from Spot 1 



Emission Line 


FWHM" 


Observed Flux'' 


Nv AA1239,1243 




2.5 ±0.9 


Si IV AA1394, 1403 




0.6 ±0.2 


Oiv] A1400 




1.3 ±0.4 


[Niv] A1483 




1.4 ±0.4 


Niv] A1487 




1.5 ±0.4 


Civ AA1548, 1551 




2.9 ±0.8 


Hen A1640 




3.5 ±0.9 


[Oi] A6300 


154 ±7 


8.9 ±0.7 


[Oi] A6364 


179 ± 14 


3.3 ± 0.4 


[Nil] A6548 


165 ± 11 


5.6 ±0.8 


Ha A6563 


253 ±4 


46.7 ± 1.2 


[Nil] A6584 


204 ±7 


13.1 ±0.9 


He I A6678 


177 ± 38 


0.8 ±0.2 


[Sii] A6717 


113 ±40 


0.9 ±0.9 


[Sii] A6731 


149 ± 48 


0.9 ± 0.4 


" in km 






in units of 10 ergs cm ^ s ^ 





In addition to interstellar extinction, the observed C IV and Si IV fluxes are further reduced by 
line absorption along the line of sight to the supernova. The absorption profiles for these lines are 
known from lUE high-dispersion spectra of SN 1987A taken at the time of outburst (Welty et al. 
1999). Since the intrinsic profiles of the lines are poorly known, correcting for this absorption is a 
somewhat uncertain endeavor. For now we just present the observed fluxes and note that under 
reasonable assumptions of the intrinsic profiles we find correction factors of 2 for these lines (Pun 
et al. 2000). 
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3. RADIATIVE SHOCK MODEL 

Figure 3 illustrates the scenario that we believe accounts for the emission from Spot 1. The 
freely expanding supernova ejecta are slowed from ^ 15, 000 km s~^ to ~ 4, 000 km by a reverse 
shock. The shocked gas drives a blast wave (forward shock), with velocity ~ 4,000 km s~^, 
into an Hii region having density nn ~ 100 cm~^ (Chevalier 8z Dwarkadas 1995; Borkowski, 
Blondin, & McCray 1997a; Lundqvist 1999). The bright spots occur where the blast wave en- 
counters protrusions of the dense (n ~ 10*^ cm"'^, Lundqvist & Fransson 1996) equatorial ring. 
As it enters a protrusion, the blast wave is slowed to a value Vg ~ vi,{nii/n)^^^ f(9), where 6 is 
the angle between the surface of the protrusion and the propagation direction of the blast wave. 
The function f{0), which accounts for shock obliquity and the pressure increase due to the re- 
flected shock, varies from /(O) ~ 2 at the head of the protrusion to /(7r/2) ~ 0.7 along its 
sides (Borkowski, Blondin, & McCray 1997b). We therefore expect a range of shock velocities, 
280 km s^ ^ ^^s ^ 800 km s~ , to be present (where 714 is the pre-shock density in 

units of 10^ cm~'^). 




Fig. 3. — Cartoon of SNR 1987A's structure in the equatorial plane. 

An important parameter for the shocks considered here is the radiative cooling time tc, defined 
as the time for the shocked gas to cool from the post-shock temperature, Tg = "ifim^Vg /IQk ~ 1.2 x 
10®(fs/300 km s~^)^ K, to Tc = 10^ K. A fit to the results from a plane-parallel steady-state (ID- 
SS) shock code (e.g.. Cox & Raymond 1985) with abundances appropriate to SN 1987A (Lundqvist 
& Fransson 1996; Russell & Dopita 1992) shows that, for 100 km s"^ < < 600 km s"\ « 
4.6nJ^(?;s/300 km s~^)3-7 yr. 

When these spectra were taken, the age of the spot, igpot ~ 3 yr, was comparable to the cooling 
time for shocks with the range of densities and velocities expected in the protrusion. Shocks with 
velocities ^ 250 n^-^^ km s"^ have > tspot and would not yet have radiated away their energy, 
while slower (radiative) shocks will have had time to cool and will convert all their thermal energy 
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into radiation. Non-radiative shocks produce far less ultraviolet and optical emission than radiative 
shocks. Therefore, although faster shocks might be present, we expect the ultraviolet and optical 
emission from Spot 1 to be dominated by shocks with speeds ^ 250 ra^'^^ km s~^. 

The radiative shock model does predict the presence of both strong UV and optical lines. The 
ultraviolet lines produced in a radiative shock come mainly from the cooling region, where the 
shocked gas's thermal energy is radiated away. The optical lines, however, are formed in the cool 
layer of dense (~ nTs/Tc, where ^ 10^ K) gas which forms behind the cooling region. This 
photoabsorption layer quickly becomes thick enough to absorb the downstream ionizing radiation 
produced in the cooling region and converts the flux into optical emission lines. 

The optical lines are formed with thermal widths characteristic of 10^ K gas. The actual 
profiles are much broader and must be due to the macroscopic motion of fluid parcels in the 
photoabsorption layer. In a steady shock, the gas in the photoionization zone will have the same 
velocity as the shock front, and therefore the optical lines would be expected to have velocity 
profiles representative of the hydrodynamics of the shock (specifically, those parts of the shock that 
have developed a radiative layer). The presence of both redshifted and blueshifted emission can be 
explained by shocks traveling both into and out of our line of sight (see Figure 3) . 

4. DISCUSSION 

It will not be an easy task to develop a quantitative model for the spectrum and line profiles 
of the hot spot. As Figure 3 illustrates, these result from a superposition of shocks having a range 
of velocities and aspect angles that depend on the unknown geometry and density structure of the 
protrusion. Our task is further complicated by the fact that radiative shocks are subject to violent 
thermal and dynamical instabilities (e.g., Innes, Giddings, & Falle 1987; Klein, McKee, &: Colella 
1994; Walder &; Folini 1996, 1998), so steady shock models may be inadequate to interpret the 
spectrum. 

Despite the complexity of the actual situation, we can gain some insight into the physical 
conditions in the spot by comparing the actual spectra with the results from the ID-SS shock code. 
For example, the code shows that the ratios of the flux of N V AA1239, 1243 to the flux of lines from 
lower ionization states (Oiv] A1400, Niv] A1487, and Civ AA1548, 1551) are sensitive functions of 
shock velocity. N V is not produced at temperatures ^ 2 x 10^ K (v^ <^ 120 km s^"*^), above which 
its emissivity rises rapidly with shock velocity. On the other hand, O iv, N IV, and C iv all form at 
lower temperatures and their emissivities are less sensitive to velocity. To reproduce the observed 
ratios with the shock code we must constrain the shock velocity to the range 120 km s~^ < Vg < 
150 km s~^. But such a velocity range conflicts with the observed optical line widths, which require 
radiative shocks with velocities up to at least ^ 250 km s~^. 

However, only 20% of the emission in the broad Ha line profile comes at high velocity 
{\v\ > 135 km s~^, measured from line center). We may then argue that only 20% of the Ha 
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emission must come from shocks faster than 135 km s~ . If this is true, the hne ratios will be 
dominated by emission from the slower shocks. Wc find that if the slower shocks {vg < 135 km s~^) 
cover about four times the surface area than the faster shocks {vg > 135 km s~^), then the net 
emission from these shocks produced ratios near those observed. 

The presence of both fast and slow radiative shocks indicates that the protrusion must have a 
slightly higher density than the average ring density. Using hydrodynamical arguments presented 
in §3 we find that preshock densities n ^ 4 x 10^ cm~^ are required to produce shocks as slow 
as 135 km s~^. A less stringent lower limit is provided by an independent argument based on 
the radiative cooling time. Preshock densities of n ^ 10^ cm~^ are required for shocks as fast as 
250 km s~^ (as indicated by the Ha line width) to cool within tg ^ S years. 

We can also use the spectroscopic data to estimate the surface area of the shock interaction. 
Most of the emission comes from shocks with velocities ~ 135 km s~^, for which the shock code gives 
an Ha surface emissivity of 0.4 71,4 ergs cm~^ s~^. Correcting the observed Ha flux for reddening 
(Scuderi et al. 1996) and assuming a distance of 50 kpc, we find that the shock must have an 
Ha emitting surface area of Ag ~ 10^^ cm^ (assuming n = 4 x 10^ cm~^). Then, assuming a 
hemispherical shape for the emitting surface, we estimate that Spot 1 should have a characteristic 
dimension Dgpot ~ (2^s/7r)^/^ ^3x 10^^ cm (or ^ 0^'04). This is smaller than HST's resolution so 
we are not surprised that the spot is spatially unresolved. 

As Table 1 indicates, different emission lines have different widths. This is expected, since the 
emissivities of the lines are sensitive functions of shock velocities, and a range of shock velocities is 
clearly in evidence. Moreover, the differences among the optical line widths may be explained by 
partial cooling in fast shocks. Emission from ions like O I and S li comes from regions downstream of 
the gas that produces Ha emission. Faster shocks take longer to develop these downstream regions. 
The fastest shocks that produce Ha will not yet have developed their [O l] and [S ll] emitting zones. 
Therefore, the [O l] and [S ll] lines are narrower than the Ha line. 

Now that several new bright spots have appeared on the circumstellar ring of SN 1987A 
(Lawrence Sz Crotts 2000; Garnavich, Kirshner, & Challis 2000; Lawrence et al. 2000), we have 
new opportunities to elucidate the development of radiative shocks. By observing the spectra of 
each of the spots, we can follow the development of shocks having different ages and velocities. The 
spots will continue to brighten and their spectra will change. The pressure that drives the shocks 
will increase, the shocks will engulf more of the protrusions, and faster shocks will develop radiative 
layers. The most powerful constraint on the hydrodynamics will be to resolve the line ratios as a 
function of Doppler velocity. For example, the wings of the N V AA1239, 1243 line profile should be 
more prominent than those of CiV AA1548, 1551. As the optical lines brighten, their wings should 
broaden. 

For these reasons, it is imperative to observe the development of the spectra of these spots on 
a regular basis. These observations will need to be performed by STIS in order to obtain UV line 
profiles as well as spatially distinguish the spectra of the different spots. 
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